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CN . ABSTRACT 

Context. This article investigates Sloan Digital Sky Survey (SDSS) star counts based on the Just-JahreiB Galactic disc model by 
analysing the TRI-dimensional modeL of thE GALaxy (TRILEGAL) and Besangon models of the Milky Way 
I Aims. We test the ability of the TRILEGALand Besangon models to reproduce the colour-magnitude diagrams of SDSS data at 

■ the north Galactic pole (NGP). We show that a Hess diagram analysis of colour-magnitude diagrams is much more powerful than 
luminosity functions in determining the Milky Way structure. 

Methods. We derive a best-fitting TRILEGAL model to simulate the NGP field in the {g - r, g) colour-magnitude diagram of SDSS 
I filters via Hess diagrams. For the Besangon model, we simulate the luminosity functions (LP) and Hess diagrams in all SDSS filters. 

OWe use a analysis and determine the median of the relative deviations in the Hess diagrams to quantify the quality of the fits 
by the TRILEGAL models and the Besangon model in comparison and compare this to that of the Just-JahreiB model. The input 
\ isochrones in the colour-absolute magnitude diagrams of the thick disc and halo are tested via the observed fiducial isochrones of 

O I globular clusters. 

I ■ Results. We find that the default parameter set lacking a thick disc component gives the best representation of the LP in TRILEGAL. 

O ' The Hess diagram reveals that a metal-poor thick disc is needed. In the Hess diagram, the median relative deviation of the TRILEGAL 

I model and the SDSS data amounts to 25 percent, whereas for the Just-JahreiB model the deviation is only 5.6 percent. The isochrone 

C/^ , analysis shows that the representation of the main sequences of (at least metal-poor) stellar populations in the ugriz filter system is 

■ reliable. In contrast, the red giant branches fail to match the observed fiducial sequences of globular clusters. The Besangon model 
shows a similar median relative deviation of 26 percent in (g - r,g). In the u band, the deviations are larger. There are significant 
off'sets between the isochrone set used in the Besangon model and the observed fiducial isochrones. 

Conclusions. In contrast to Hess diagrams, luminosity functions are insensitive to the detailed structure of the Milky Way components 
due to the extended spatial distribution along the line of sight. The flexibility of the input parameters in the TRILEGAL model is 
insufficient to successfully reproduce star count distributions in colour-magnitude diagrams. In the Besangon model, an improvement 
of the isochrones in the ugriz filter system is needed in order to improve the predicted star counts. 
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o: 

■ 1 - Introduction ing photometric data in five bands of 2.6 x 10^ stars and spectra 

■ I 1 J I . of more than 50 0,000 stars (for the most recent eighth data re- 

• • Since the seminal work of [ Bahcall & Soneira| ( |1980a|H 11984D, lease (DR8), see lAihara et al .112011 . The SDSS has become one 

^ the method of simulatmg star counts based on analytic Milky ^^^^ ^i^^l^ ^.^^ databases for Galactic and extragalactic 
[Xj Way models has evolved considerably. The structural and evolu- astronomy in recent years. The availability of SDSS data enables 
tionary parameters of the Milky Way components have been re- ^lake significant progress in constraining Galactic model 
5t ^ fined in multi-population scenarios. However, a "concordance" parameters 
model including the detailed shape and parameters of the den- 
sity profiles, star formation history (SFH), age-velocity disper- 
sion relation (AVR), etc., has not yet been reached. Analytic models are a powerful tool to constrain either evo- 
On the other hand, there is a rapidly increasing set of obser- lutionary scenarios or assumptions about Galactic structure by 
vations consisting of more accurate and much deeper photome- comparing model predictions with a large variety of observa- 
try to constrain the parameters of the models of each population, tional constraints such as star counts and kinematics. Models can 
The observations include wider and deeper sky coverage and ad- be used to create probability distributions in the parameter space 
ditional filter systems providing results that may require modifi- of observable or mock catalogues via Monte Carlo simulations, 
cations of the existing detail ed Milky Way models. The Sloan Realistic parameters and details for further improvement can be 
Digital Sky Survey (SDSS, [York et al.| |2000|; [Stoughton et^ap obtained by comparing model predictions with suitable observa- 
|2002|) covers more than one third of the celestial sphere provid- tions, such as multi-directional photometric star counts. Within 
these parameters and functions, the density profiles, SFHs, and 
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To simulate and predict star counts in every direction of the 
sky in the bands of various photometric systems, several auto- 
mate d tools are av a ilable on the web, e.g., the "Besan9on" model 
Q by iRobin et al.l (l2003h and the TR I-dimensional modeL of 
thE GALaxy ("TR ILEGAL") model0 (iGirardi & Bertellill 19981: 
iGirardi et al.l[2Q02[l2QQ5h . which are able to deal with isochrone 
construction, synthetic photometry, and simulations of t he stel- 
lar pop ulations of the main Galactic components. Sharma et al.l 
(I2OIII) presented Galaxia, a new tool which combines the ad- 
vantages of the Besan^on and the TRILEGAL model. The code 
generates very efficiently synthetic surveys and allows one to 
choose between a much larger variety of input models. The first 
version of Galaxia is expected to be presented online soonB 

The "Besan^on model" (iRobin et al.ll2003l) is presently the 
most elaborate and refined publicly available tool for predict- 
ing star counts. In this model, the Milky Way is divided into four 
components (thin disc, thick disc, spheroid, bulge), which are de- 
scribed by their spatial distribution, SFHs, initial mass functions 
(IMFs), sets of evolutionary tracks, kinematics, and metallicity 
characteristics, and include giants and a white dwarf population. 
From Monte Carlo simulations, mock catalogues including ob- 
servable parameters as well as theoretical ones are obtained. 

The TRILEGAL model is a Monte Carlo simulation to pre- 
dict the probability distributions of stars across small sky areas 
and for special passbands systems, is still work in progress. It 
is a population synthesi s code to simulate photometry of any 
field in the Milky Way (IGirardi et al .1 120051). It was developed 
to model population synthesis star counts of the Milky Way. It 
allows its users to create a pseudo stellar catalogue including po- 
sitions, multi-colour photometry, and other physical parameters 
according to user-defined structural parameters. 

On the basis of SDSS photometry from Data Release 6 
(DR6, lAdelman-McCarthv et al.|[2008h . lJuric et al.1 (l2008h de- 
rived a three-dimensional (3D) stellar distribution of the Milky 
Way by applying universal photometric parallaxes. They then 
fitted two discs, each with exponential profiles in the radial and 
vertical directions. They found vertical scale heights of 300 pc 
and 900 pc for the thin and thick disc, respectively. The lo- 
cal density normalisation of the discs and basic parameters of 
the stellar halo were also determined. Subtracting the resulting, 
smooth Milky Way model, overdensities (and underdensities) 
were revealed that could be investigated in more detail. 

We have developed a new local Milky Way model based 
on local s tellar k inematics and star counts in the solar neigh- 
bourhood (lJust & J ahreiB 2010, hereafter "Just-JahreiB model"). 
Compared to the Besangon and TRILEGAL models, the main 
advantages of our approach are that vertical density profiles are 
fully consistent with the SFH and AVR in the gravitational field 
and that we use a well-calibrated empirical main- sequence (MS). 
In the current state of the Just-JahreiB model, only MS stars 
are used instead of including the full Hertzsprung-Russell dia- 
gram for each stellar sub-population. The comparison of the ba- 
sic Just-JahreiB disc model with SDSS star counts in the NGP 
field with b > 80° en abled us to select a best-fitting SFH and 
AVR (I Just et al.l l201ll) . Predicted star-number densities fit the 
data with a median relative deviation of 5.6% across the Hess di- 
agram. We use the Just-JahreiB model to judge the ability of the 
TRILEGAL and the Besangon model to reproduce SDSS data at 
the NGR 



^ http://model.obs-besancon.fr/ 

^ http : //stev. oapd. inaf . it/cgi-bin/trilegal 

^ http://galaxia.sourceforge.net/ 



In this article, we present a test of the ability of the 
TRILEGAL code to reproduce the diff'erent Milky Way mod- 
els mentioned above, since the TRILEGAL model is presently 
the only code that allows an interactive analysis of the im- 
pact of model parameter variations on star count predictions. 
Additionally, we try to determine a best-fit model by a system- 
atic input-parameter optimisation in the TRILEGAL model. The 
Besangon model has no free parameters and we perform a multi- 
colour analysis over all five filters in the ugriz system. We restrict 
the analysis to the SDSS star-count data of the NGP field, since 
these are very high quality data and they give a direct measure 
of the controversial vertical structures of thin and thick discs. 

The output catalogues of the TRILEGAL and Besan^on 
models with photometric parameters can be used to analyse the 
contributions of diff'erent populations to the observed star counts 
and to quantify systematic deviations in order to identify the 
weaknesses that need to be resolved. 

Sect. 2 introduces the data used in this work. In Sect. 3, we 
describe the details of the Just-JahreiB, TRILEGAL, Besangon, 
and the Juric models used in our analysis. In Sect. 4, we present 
the results of the TRILEGAL model concerning the isochrones, 
luminosity functions, and Hess diagrams. In Sect. 5, we present 
the comparisons between the Besangon model and the data us- 
ing similar methods. In Sect. 6, we summarise the results of our 
study. 

2. Data 

The SDSS has become one of the most widely used databases for 
Galactic and extragalactic astronomy in recent years. The eff'ec- 
tive exposure time for imaging is approximately 54 seconds with 
a limiting magnitude of r ~ 22.6 mag. The g and r band point- 
source imaging of the SDSS ranges from 14 mag to 22.5 mag. 
We selected stars from the Catalog Arch ive Server (CAS) of th e 
seventh data release of the SDSS (DRT. lAbazaiian et al.ll2QQ9h . 
The SDSS table Star, which we utilise, contains the photometric 
parameters (no redshifts nor spectroscopic parameters) of all pri- 
mary point-like objects classified as stars from Photo Primary . 
We only consider the NGP field (Galactic latitude /? > 80°) in 
this publication. A series of lower latitude fields will be anal- 
ysed in the future. 

We use de-reddened magnitudes from the SDSS point- 
spread function phot ometry b ased on the values from the ex- 
tinction map of >Schlegel et al.l ([199 81) . In the following, all mag- 
nitudes and colours are de-reddened. We select all stars in DR7 
within a colour range of (g-r) = [-0.2, 1 .2] and magnitude lim- 
its of g = [14, 20.5] with magnitude errors below 0.2 mag in g 
and in r. In these magnitude and colour ranges, the error limits 
remove 0.88% of the sample in the NGP field. A total number of 
N = 276 180 stars remain in the final sample. 

For our analysis, we produce Hess diagrams of the stel- 
lar number-density distribution in colour-magnitude bins nor- 
malised to 1 deg^ sky area in the colour-magnitude diagram 
(CMD) for a given direction. One obvious deficiency are the un- 
known stellar distances, hence the vertical axis of the CMDs nec- 
essarily shows apparent magnitudes instead of absolute magni- 
tudes. The combined distribution of (apparent) stellar luminosity 
and colour shows the spatial structure and the population prop- 
erties in the chosen direction, which can then be compared with 
model-generated Hess diagrams. 

In the top panel of Fig. [T] the Hess diagram in (g - r, g) is 
shown. Number densities per square degree, 1 mag in g and 0.1 
mag in (g-r), are colour-coded on a logarithmic scale. The bin 
size is (0.05, 0.5) in (g - r, g) with additional smoothing in steps 
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Table 1. Properties of five fiducial globular clusters observed in 
the SDSS filter system. 
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Fig. 1. Top: Hess diagram of the NGP stellar sample taken from 
SDSS DR7 and covering 313.36 deg^. The x-axis is the de- 
reddened colour index (g - r) and the y-axis is the de-reddened 
apparent magnitude g. Number densities per square degree, 1 
mag in g and 0.1 mag in - r), are colour-coded on a log- 
arithmic scale. The stellar number density has been smoothed. 
Bottom: De-reddened g magnitude distribution of stars in the 
same colour range but with an extended magnitude range. The 
bins of this luminosity function are the logarithmic star numbers 
per 0.5 mag. The assumed Poissonian error bars are also indi- 
cated. (A colour version is available on-line.) 



of 0.01 mag. From the upper right (purple/blue regions) to the 
lower left (red triangular regio n) the thin disc, thick disc, and 
halo dominate, respectively (see I Just et al.ll201 iT) . In the top left 
of the figure, no stars are visible. In the bottom left, white dwarfs, 
blue horizonal branch stars, and blue stragglers of the thick disc 
and halo may appear at (g - r) < 0.2, but the majority of objects 
are mis-identified point-like extragalactic sources. 

The luminosity function is obtained by projecting the Hess 
diagram along the (g - r)-axis. In the lower panel of Fig.[Tl the 
SDSS data are shown for the complete range in the g band in or- 
der to illustrate the total numbers of stars and the cutoff's at bright 
and faint magnitudes. Poissonian error bars were adopted. The 
incompleteness at g > 22 mag is obvious and a significant con- 
fusio n with m is-identified extragalactic sources is likely to occur 
here (I York et al. 2000). To be on the safe side, we choose for our 
analysis g = 20.5 mag as the faint limit. One significant jump 
appears at g = 14 mag, which is caused by the incompleteness 



ID 


(m - M)o 


E(B - V) 


[Fe/H] 


M71 


13.02 


0.25 


-0.73 


M5 


14.37 


0.03 


-1.27 


M13 


14.42 


0.02 


-1.54 


NGC4147 


16.42 


0.02 


-1.83 


M92 


14.58 


0.02 


-2.28 



Notes. Col. 1 lists the names of the globular clusters. Cols. 2 and 3 con - 
tain distance modulus and reddening E( B - V) from | An et all (l2008h . 
Col. 4 lists the metallicity of each cluster. iHarri si (1 1996b metallicitv val- 
ues are adopted here to ensure that we have a consistent set. These 
clusters are used to enable us to compare with the TRILEGAL and 
Besan9on models. 



due to the large magnitude errors. The bright limit at 14 mag is 
set by the saturation of the stars in the SDSS observations. 

Apart from the DR7 catalogue data, we select the fidu- 
cial is ochrones of five globular clusters determined by lAn et al.l 
(I2OO8I) to compare the photometric properties of the thick disc 
and halo populations in the TRILEGAL and Besan^on models. 
Using the extinct ion paramete rs and distance moduli of these 
globular clusters (lHarrislll996b , the de-reddened absolute mag- 
nitudes and colour indices are determined. The parameters of the 
globular clusters are listed in Table [T] These clusters are selected 
to cover the metallicity range of the thick disc and halo adopted 
in the TRILEGAL and the Besan^on models. 



3. Models 

We describe the basic properties of the four models compared 
with the SDSS data. 



3.1. Just-JahreiB model 



lJust & JahreiBi (1201 Ol) presented a local disc model where the 
thin disc is composed of a continuous series of isothermal sub- 
populations characterised by their SFH and AVR. A distinct 
isothermal thick disc component can be added. The vertical den- 
sity profiles Ps,j(z) of all sub-populations j are determined self- 
consistently in the gravitational potential 0(z) of the thin disc, 
thick disc, gas component, and dark matter halo. The density 
profiles of the MS stars are calculated by adding the profiles 
of the sub-populations, accounting for their respective MS life- 
times. In a similar way, the thicknesses and the velocity distribu- 
tion functions of MS stars are derived as a function of lifetime. 

For the thin disc of the Just-JahreiB model, a factor of 
= 1.53 with geff = 0.654 needs to be applied in order 
to convert the present-day stellar density to the density of ever 
formed stars. For the thick disc and halo, the conversion factor 
is 1.8.^ 

In I Just & JahreiBi (l2010l) . we showed that the vertical profile 
of an isothermal thick disc can be well-fitted by a sech^^(z/at^t) 
function. The power-law index at and the scale height ht are con- 
nected via the velocity dispersion. 

The SFH and the AVR are the main input functions that are 
to be determined as a pair of smooth continuous functions. For 
each given SFH, the AVR is optimised to fit the velocity dis- 
tribution functions of MS stars in V-band magnitude bins in 
the solar neighbourhood. A metal enrichment law (i.e., the age- 
metallicity relation (AMR)) is included that reproduces the local 
G-dwarf metallicity distribution. Additionally, a consistent IMF 
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is determined, that re produces the local luminosity function. In 
lJust & JahreiBi (|201Q|) . four different solutions with similar min- 
imum but very different SFHs were discussed. 



In Just et al 



(1201 Ih . we combined the models of 

I Just & JahreiBi (|201Q|) with the empirically determined photo- 
metric properties of the local MS in the ugriz filter system. 
By adopting a best-fitting procedure to determine the local 
normalisations of thin disc, thick disc, and stellar halo as a 
function of colour, we compared the star count predictions of the 
different models for the NGP field. In this procedure, the thick 
disc parameters at and ht were optimized and we added a simple 
stellar halo with a flattened power-law-density distribution. 

To measure the quality of the model and fit the data, we here 
define slx^ value for the Hess diagram. Each Hess diagram con- 
tains mxn boxes of magnitudes and colours. The relative differ- 
ence (data - model)/model of each Hess diagram is a matrix Dij 
with the indices / for colour and j for magnitude. 

The total is used to measure of the modelling quality, 
which is our standard way of evaluatin g the goodness of a model. 
We apply the same definition of as lJust et aP (l201ll) to mea- 
sure the goodness of the TRILEGAL and Besan9on models. We 
also define the median deviation Amed = median(|Z)/^|) of the 
relative differences Dij for each Hess diagram to quantify the 
typical difference between model and data. 

In lJust et al.1 (l201ll) it turned out that the lJust & Jahreil 
(I2QTqI) "model A" is clearly superior to the other models. We 
therefore select "model A" as the fiducial Just-JahreiB model. Its 
basic density profiles and parameters are listed in the last column 
of Table [3] for reference. 



3.2. TRILEGAL 

The main goal of the TRILEGAL model is the simulation of 
Milky Way star counts in several p hotometri c system s using sim- 
ple profiles. For this purpose, Gi rardi et"aD (l2002h developed a 
detailed method to convert theoretical stellar spectral libraries 
to photometrically measurable luminosities. These authors can 
produce these transformations routinely for any ne w photometric 
system , such as the five-colour filters of the SDSS (iGirardi et al.l 
|2004|) . In addition, TRILEGAL has been used to predict stellar 
photometry in different systems. 

Being able to simulate both very deep and very shallow sur- 
veys is one of the most difficult goals of TRILEGAL. For shal- 
lower surveys such as Hipparcos data, MS, red giant branch, and 
red clump stars with partly very low number statistics need to 
be generated. On the other hand, for more sensitive surveys an 
extended lower-luminosity MS, which 

The measurements of stellar age (r) and metallicity (Z) pro- 
vide parameters such as bolometric magnitude Mtoi, effective 
temperature Teff , surface gravity g, core mass, and surface chem- 
ical composition. 

The IMF is a crucial ingredient in Galactic modelling, be- 
cause it determines the relative number of stars with different 
photometric properties and lifetimes. For comparison and con- 
venience, the default sha pe of the IMF is a log-normal function 
in TRILEGAL, following Ichabrid (l200lh . 

The SFH describes the age distribution of the stars. For old 
components, the shape of the SFH is quite poorly constrained, 
because the observed star counts are quite insensitive to it. For 
the thin disc, the SFH has not yet been well-determined (e.g., 
lAumer & BinnevI (l2009l) : lJust &Jah reiB (2010)). The interface 
of TRILEGAL allows only the choice between a constant star 
formation rate (SFR) up to an age of 1 1 Gyr and a two-step SFH 



with a L5 times enhancement in the SFR between the ages of 1 
and 4 Gyr. 

The AMR is also an important ingredient of the Galactic 
model, because it can h elp modify the stellar CM Ds. For the thin 
disc, the abundances of iRocha-Pinto et aP (l2000l) are adopted in 
TRILEGAL. 

There are five Galactic components in TRILEGAL: thin disc, 
thick disc, halo, bulge, and an extinction layer in the disc. 

The mass density of the thin disc decreases exponentially 
with Galactocentric radius R projected onto the plane of the disc. 
The Galactic thin disc is radially exponential and locally not 
isothermal, because the scale height is a function of stellar age 
and the thin disc is not isolated. In TRILEGAL, one can choose 
between an exponential and a sech^ profile for the vertical densi- 
ties of both the thin and the thick disc. The typical height of the 
disc sub-populations depends on stellar age. The stellar density 
in the solar neighbourhood pd(0) is the mass density of all stars 
ever born there. 

The density of the thick disc can be described as either a 
double exponential or vertically with a sech^ form similar to 
that of the thin disc. However ht is not a function of stellar age 
because the thick disc is predominantly an old population. The 
SFH is constant over an age range of 11-12 Gyr with a metal- 
licity [Fe/H] = -0.67 ±0.1 dex and an a-enhancement of [a/Fe] 
-0.3. 

The mass density Ph of the halo follows a de Vaucouleurs r^^"^ 
law (de Vauc ouleursI 1 94 8l) and allows for a spheroidal flattening 
(iGilmore.J984 . The density profile of the stellar halo can be de- 
scribed by the parani eters effective radius, oblateness, and local 
density (lYounglll976l) . The SFH is constant over an age range of 
12 — 13 Gyr and there are two choices for the metallicity distri- 
bution. We fix the metallicity to be [Fe/H] = -1.6 ± 1 dex with a 
corresponding of-enhancement of [a/Fe] = 0.3. 

The bulge has not yet been calibrated well in TRILEGAL be- 
cause of a lack of observational constraints. Since the bulge does 
not contribute to high Galactic-latitude data, we do not discuss 
its parameters any further. 

We do not use the dust extinction model of TRILEGAL but 
instead we compare star counts with de-reddened observational 
data as described in Sect. 2. 

Each simulated photometric catalogue in TRILEGAL has a 
given direction that is a conical area with a maximum sky cover- 
age of 10 deg^. The centre of the area can be specified in either 
Galactic coordinates (/, b) or equatorial coordinates (a, S). The 
resolution in magnitudes. Am, in TRILEGAL's calculations can 
be modified by the user. The default value is Am = 0.1 mag and 
can be reduced to a minimum of 0.05 mag. Any data point with 
a resolution less than Am will not be presented. The output cata- 
logues of TRILEGAL are based on a random number generator, 
which allows one to investigate the effect of noise. 

Since the maximum sky area is much smaller than the ob- 
served field (313.36 deg^ for the NGP field), we generate cat- 
alogues for 37 different sight lines (Table [2]) to cover the NGP 
field. We add up all the star counts and then normalise to star 
counts per deg^. 

3.3. Besangon model 

iRobin et aP (l2003l) developed the population synthesis approach 
to simulate the structure and evolution of the Milky Way. The 
Besangon model is based on different well-calibrated data sets 
and the structural parameters are optimised to reproduce a set of 
selected fields on the sky. The construction of each component 
is based on a set of evolutionary tracks and assumptions about 
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stellar number-density distributions, which are both constrained 
by either dynamical considerations or empirical data and guided 
by a scenario of a continuous formation and evolution of stellar 
populations. 

The Besangon model deals with the thin disc is com- 
posed of seven isothermal sub-populations with different age 
ranges, kinematics parameters (W velocity dispersion ctw), scale 
heights, and local densities. The thick disc is considered as a 
single population with an age of 11 Gyr and a local density 
value (1.34 x lO'^Mopc"^). The halo is modelled by a sin- 
gle population of older age (14 Gyr) and lower local density 
(2.2 X 10"^Mopc"^ including white dwarfs). This prescr iption 
is similar to our Just-JahreiB model (I Just & JahreiBll20l"Qb . The 
main differences are that 1) in the Just-JahreiB model, we adopt a 
continuous series of isothermal sub-populations in the thin disc; 
2) the density profile and local density of each component are 
determined by solving iteratively Poisson and Jeans equations 
that are consistently based on the local kinematic data and the 
SDSS photometric data; and 3) we use a mean empirically cali- 
brated MS instead of a full Hertzsprung-Russell diagram (HRD) 
based on stellar evolution. 

Since we limit our targets to the NGP field with a Galactic 
latitude b > 80°, stars belonging to the Galactic bulge can be 
ignored. The stellar populations in the NGP field can be divided 
into three distinct components: the thin disc, the thick disc, and 
stellar halo. For d etails of the propert ies of the Galactic compo- 
nents, we refer to iRobin et al.1 (l2003h . With the NGP field, we 
can test the vertical structure of the Galactic components and 
the calibration of the underlying isochrones in the SDSS filter 
system. 

The Besan9on model outputs u^g'/i'z' photometry. To ob- 
tain SDSS magnitudes in g, r, /, and z and the colours u - g, 
g - r, r - i and / - z from the Besan^on model, we transform 
the photometric system from the u*, g\ r\ i' and z' filter set 
using a two-step method . In the first step , we transform from 
the special u* filter to u' (IClem et al.ll2008l) . In the second step, 
we transform from the u\ g\ r\ i\ a nd z' filters t o u, g, r, /, 
and z followin g the t ransformations in iTucker et"aD (l2006l) and 
lJust & Jahreil (I2008h . 



3.4. June model 



In I Juric et al ] (I2008h . a map of the 3D number density distribu- 
tion in the Milky Way derived by photometric parallaxes that 
were based on star count data of SDSS DR6 is presented. The 
stellar distances in the sample were determined by the colour- 
absolute magnitude relation of MS stars. Most stars are at high 
latitudes and cover a distance range from 100 pc to 20 kpc and 
6500 deg^ on the sky in the sample. 

The density distribution is fitted using two double- 
exponential discs and the sample implies the existence of an 
oblate halo. The estimated errors in the model parameters are 
smaller than ~ 20% and the errors in the scale heights of discs 
are smaller than ~ 10%. By adjusting the fraction of binaries 
and the bias correction, the scale heights of the discs are found 
to be /zd = 300 pc and = 900 pc, and the local calibration 
is pt(0)/pd(0) = 0.12. The best- fit of the spheroid leads to an 
oblateness of 0.64, a power-law radial profile of ph oc r"^-^, and a 
local calibration of Ph(0)/Pd(0) = 0.005. After subtraction of the 
smooth background distribution, there are still overdensity areas, 
in particular areas that a re explained as tid al streams including 
the ''Vi rgo overdensity" ([Puffau et al.ll2006l) and the Monoceros 



Table 2. The central Galactic coordinates of 37 selected fields 
simulated using TRILEGAL in order to be compared with SDSS 
observations. 

b (deg) / (deg) 
90 

87 0, 60, 120, 180, 240, 300 

84 0, 30, 60, 90, 120, 150, 180, 210, 240, 270, 300, 330 

81 0, 20, 40, 60, 80, 100, 120, 140, 160, 180, 200, 220, 

240, 260, 280, 300, 320, 340 



The solar offset Zq is (25 ± 5) pc or (24 ± 5) pc at the bright 
and faint ends, respectively, which is consistent with the value of 
24.2 pc adopted by TRILEGAL. 

4. Simulations and results of TRILEGAL 

The model TRILEGAL can be accessed via a web-based inter- 
face. Users can simulate the stellar populations in a special vol- 
ume defined by the luminosity range and direction on the sky us- 
ing different parameter sets. Because the parameter sets and pro- 
file shapes are relatively simple and the web-based programme 
can return the results in real time, we can compare the output 
interactively with the data. 

In general, it is difficult to find the "best" input parame- 
ters for the components by using the web interface interactively, 
since the transitions of thin to thick disc and thick disc to halo in 
the CMD are not known a priori. We cannot search indepen- 
dently for the best-fitting parameters in each component. We 
therefore first test both the parameters of the default model and 
those that reproduce best the other three Milky Way models de- 
scribed in Sect. O 

We fix some parameters and input functions in the 
TRILEGAL model in order to facilitate the comparison with ob- 
servational data. We fix the IMF to a Chabrier log-normal func- 
tion and use the default binary distribution with a binary fraction 
of 0.3 and a mass-ratio range of 0.7-1 . The solar position is at the 
default value of (Rq, zq) = (8700, 24.2) pc and we fix the SFH 
for the thin disc to be constant. For the other input parameters 
in the TRILEGAL model, we selected five sets of parameters as 
given in Table[3] The "set 1" parameters in this table are from the 
default input of TRILEGAL in the web interface, "set 2 " is best 
adjusted to the Just-Jahr eiB model (Just & JahreiBll201Qb . "set 3" 
to the Besangon niodel (iRobin et alll2003h . "set 4" to the Juric 
model (I Juric et al. I I2008I) . and "set 5" is optimised (as described 
in Sect. 112]). 

Because in the TRILEGAL model, we have to choose be- 
tween either an exponential or a sech^ shape of the disc profiles 



Ptiz) = pt(0)exp(--), 
Ptiz) = A(0)sech2(-^), 



(1) 
(2) 



tne Vi rgo overaensitv ( | L )i 
stream (lYannv et al.ll2003ir 



we determine for all combinations the minimum x'^ values and 
then select the profiles with the lowest value. 

At colours bluer than g - r = 0.2 mag, the SDSS data are 
unreliable owing to contamination by extragalactic objects and 
all selected models are unable to generate sufficient stars in this 
region. Therefore, the mean;^^ value (x^) of the reduced colour 
range 0.2 < g - r < 1.2 is treated as the standard way of de- 
termining the quality of the model. For completeness, the values 
(X^y for the full colour range are also listed in Table [3l 
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Table 3. The sets of input parameters in TRILEGAL and the;^^ values of the Hess diagram comparisons. 



V^U-dllLiLj 


Unit 


Set 1 


Set 2 


Set 3 


Set 4 


Set 5 


-DCSxXll^VJll IIIULICI 


Tn Qt- TjiViTf 1 R moHpl 

J U-S L J dill CIJJ IIIULICI 






IvCidUlL 


Til cf _ Tci rTrP»1 K ITTTMlf 

J UsL-J alllClJj llipUL 


JjCsall^Ull llipUL 


Til fir" iTiT^iif 
J UllC llipUL 


iiT^f 1 tni v cif 1 r\Y\ 
wpLlllllZ/ClLlUll 






\X / 




/ J. / u 


98Q 07 




94^ f\A 


908 


S4Q 
D'-ty .yo 




vt / 




151.98 


214.65 


207.60 


137.88 


104.19 


408.15 


4.31 


A J 




VJ. J)VJJ) J 






4fS^Q 


9S7rS 


0.2724 


OS^S^ 






95.00 


60.64 


13.61 


300.00 


90.00 






tn 


Gyr 


4.4 


0.7118 


0201 6 




4.4 






(y 




1.6666 


0.7576 


0.5524 





1.3 






^d,iTiax 


DC 


767 


506 


443 


300 


459 




670 


Pd(0) 






0.0297 


0.030 


0.029 




0.0394 


0.037 






59.00 


49.36 


25.20 


17.40 


50.00 




29.4 


form 




sech^ 


sech^ 


exp 


exp 


sech^ 


Einasto 




h 


pc 




793 


910 


900 


750 




800 


Pt(0) 


Mopc"^ 




0.0044 


0.0015 


0.0034 


0.0014 


0.00134 


0.0022 




Mopc-2 




7.08 


2.82 


6.12 


4.2 




5.3 


form 






exp 


exp 


exp 


sech^ 


cored exp 


sech^-^^ 


^effMh 


pc/- 


2800 


11994 


39436 


18304 


22000 


-2.44 


-3.0 






0.65 


0.65 


0.71 


0.59 


0.61 


0.76 


0.7 


Ph(0) 


Mopc-3 


0.00015 


0.000158 


0.00009 


0.000145 


0.0001 


0.00000932 


0.00015 



Notes. Cols. 1 and 2 list the parameters of each Galactic component in TRILEGAL and their units. C ols. 3 to 7 list the adopted values of the five inp ut sets for TRILEG AL. Set 1 is the default input 
of the web-based TRILEGAL interface. Set 2, 3, and 4 are from TRILEGAL fits with three models: lJust & Jahre iB (2010), Robin et al. (2Q03|), and lJuric et aP (l2008h . Set 5 (the last column) hsts 
our best input set determined by the optimisation based on the results of Sets 1 to 4. We provide two mean;^^ values and one median absolute relative diff'erence of the Hess diagram comparisons of 
the five input sets: (x^Y is with all stars covering the full colour range in the Hess diagram, (x^) is only stars with 0.2 < g - r < 1.2, A^ed gives each median value of the absolute relative diff'erence 
between the data and the model covering only for stars with 0.2 < g - r < 1.2 in the Hess diagram. The scale-height parameters of the thin disc and the local thin disc's density pd(0) as well as 
the surface density Sd(0) are listed, hi is the scale height of the thick disc. pt(0) and 2t(0) are the locally calibrated volume and surface densities of the thick disc. The reff is the eff'ective radius of 
the halo along the major axis with an oblateness ^h- The ofh is the index in the power law of the halo profile. One value of rgff and ofh is adopted based on the form of the halo profile. The ^hCO) is 
the local volume density of the halo. Cols. 8 and 9 list the corresponding parameters of the Besangon and the Just- JahreiB models for reference. 
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Fig. 2. The vertical density profiles of thick disc (top panel) and 
halo (bottom panel) of parameter set 1 (black, thick solid line), 
set 2 (red, solid line), set 3 (blue, dashed line), set 4 (green, dot- 
ted line), and set 5 (cyan, dashed-dotted line). There is no thick 
disc in set 1 . The z ranges of the panels are diff'erent. The vertical 
axis pt and ph are the densities of the thick disc and the halo in 
units of Mqpc"^. (A colour version is available on-line.) 



Fig. [2] compares the vertical density profiles of the thick disc 
and halo for the diff'erent sets. There are significant diff'erences 
between the models, even in the distance regime where the thick 
disc component dominates the star counts in the Hess diagram. 
We discuss this point later. 

To separate the impact of the synthetic colours and lumi- 
nosities of the underlying stellar populations from the struc- 
tural parameters of the Milky Way model, we compare the 
colour-absolute magnitude diagrams provided by TRILEGAL 
with the observed fiducial isochrones of globular clusters with 
corresponding metallicities. This comparison is only useful for 
both thick disc and halo, which are essentially simple stel- 
lar populations with small ranges in age and metallicity. The 
colour-absolute magnitude diagrams of Set 2 are plotted in 
Fig. [3 The cluster M71 ([Fe/H] = -0.7 3) is u sed to provide 
fiducial isochrones of the thick disc (see lAn et al. 2008). The 
three clusters M71, Ml 3, and M92 are used to provide fiducial 
isochrones of the ha lo whose mean metallicity is given as -1.6 
(iGirardi et al.l l2005h . The consistency between the model and 
the fiducial isochrones for the location of the MS and the turn- 
off' implies that the stellar population synthesis of TRILEGAL 
is satisfactory for predicting the dominant contribution to the 
star counts. However, there are some systematic inconsisten- 




0.4 0.6 

g-r (mag) 




Fig. 3. The fiducial isochrones of globular clusters from lAn et al.l 
('2008') are overplotted in the Hess diagrams of absolute mag- 
nitude. The upper panel shows the thick disc simulated by 
TRILEGAL using parameter set 2. The fiducial isochrone of the 
globular clusters M71 is overplotted as a white line. The lower 
panel shows the corresponding halo populations. From left to 
right, white lines represent the fiducial isochrones of the three 
globular clusters M71, M13, and M92 (Table [B covering the 
halo's metallicity range. (A colour version is available on-line.) 



cies concerning the location of the giant branches of the stellar 
isochrones. 



4.1. Default input (set 1) 

The default input of TRILEGAL is adopted directly as parameter 
set 1. The parameters are listed in Col. 3 of Table [3l The default 
set 1 has no thick disc. For the vertical structure pointing to the 
NGP, set 1 simulates star counts using only thin disc and halo. 

The top row of Fig. |4] shows the luminosity functions in the 
g and r bands of the SDSS data and of the models in the colour 
range 0.2 < g - r < 1.2. The default model (set 1) fits the data 
very well in the selected magnitude range 14-20.5 mag. 

In the CMD, the picture is diff'erent. We construct a Hess di- 
agram that is similar to the observed Hess diagram in Fig.[TJ The 
top row of Fig. \5\ shows the analysis of the Hess diagram. The 
left panel shows the simulated Hess diagram, which we compare 
directly to the observed Hess diagram. The middle panel shows 
the median absolute relative diff'erences with a colour-coded 
range of ±0.5. The right panel gives the unreduced distribu- 
tion on a logarithmic scale. From the relative diff'erences, some 
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Fig. 4. The magnitude distributions of the SDSS observations and the TRILEGAL simulations generated with parameter sets 1 to 4 
(top to bottom). The x- and y-axes are similar to the bottom panel of Fig. [T] The error bars (blue) are from SDSS DR7 data with the 
colour limit of 0.2 <{g-r) < 1.2. The black solid line represents the model star count predictions of TRILEGAL. The dotted (red), 
dashed (cyan), and dot-dashed (green) lines stand for the thin disc, thick disc, and halo of the TRILEGAL model, respectively. (A 
colour version is available on-line.) 
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-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Fig. 5. Hess diagrams (first column), relative diff'erences (second column) between the Hess diagrams of the models and the SDSS 
observations, and the diagrams (last column). From top to bottom, we show the comparisons based on the parameter sets 1 
to 4, respectively. The x- and y-axes correspond to the colour index -0.2 < g - r < 1.2 and de-reddened apparent magnitude 
14 < g < 20.5. The left column shows the stellar density using a colour coding from to 100 stars per square degree per 1 mag and 
0.1 colour mag on a logarithmic scale. The middle panels are the relative diff'erences between data and models with a colour coding 
from -50% to +50%. The right panels are the In;^^ distributions with a colour coding f rom -2 to +2. Note that the colour ranges in 
the middle and right panels cover much larger ranges than in Fig. 5 o f lJust et al.l (I201lh . (A colour version is available on-line.) 



clear features can be identified (middle panel of the top row of 
Fig.O: 

1. The large white "island" in the middle of the diagram is a 
consequence of the lack of a thick disc in the model. 

2. The missing stars at the top-right corner of g - r ~ 1.1 and 
^ ~ 15 mag are a hint of an underestimation of the local thin- 
disc density caused by the strongly flattened sech^ profile 
close to the Galactic plane. 

3. The underestimated number counts at the faint blue end of 
0.2 < g - r < 0.4 and g > 19.5 mag show that the halo 
density profile falls off' too steeply. 

4. The horizontally adjacent black-and-white coloured areas in 
the colour range 0.15 < g-r < 0.45 at g ~ 17 mag show that 



the (thick disc) population is replaced in the model by a more 
metal-poor (halo) population with a bluer F-turnoff' region. 
This is additional strong evidence that a thick disc compo- 
nent of intermediate metallicity is missing in the model. 
5. At the faint end (0.15 < g - r < 0.25 and 15 < g < 19), the 
stars missing in blue colours (halo) are compensated for by 
too many (or slightly too blue) M dwarfs of the thin disc in 
the luminosity function. 

The value of the mean reduced (x^) = 15 1 .98 is significantly 
larger than that of the Just-JahreiB model (= 4.31, Table O. The 
median relative deviation of the data from the model in the Hess 
diagram Amed is 0.3035, which is 5.4 times larger than for the 
Just-JahreiB model. 
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A comparison of Figs. |4] and \5\ indicates that the analysis 
of the luminosity function integrated over a large colour range 
alone may result in misleading conclusions about the underlying 
stellar populations. 



4.2. Just-JahreiB input (set 2) 

For parameter set 2, we derive the best representation of the input 
function of the Just-JahreiB model with the TRILEGAL parame- 
ters. With this "model of a model", we investigate the capability 
of TRILEGAL to reproduce a Milky Way model suggested by 
other authors. 



4.2.1. Input parameters 

For parameter set 2, which is based on the Just-JahreiB model, 
we derived a best-fit TRILEGAL profile of the thick disc to the 
sech^^ profile. Both the exponential and the sech^ profiles are 
fitted with a standard routine and the unreduced values are 
used to select the best fit. The exponential form results in = 
4.0 X 10"^, which is a better fit than the sech^ profile with;^^ = 

6.6 X 10"^. We therefore use the exponential profile. The scale 
height of the thick disc is 793 ± 6 pc and the local density is 
0.0029 + 0.000015 Mopc'l The derived local surface density of 
the thick disc is XtCO) = 4.63 Mopc"^. 

In TRILEGAL, the scale height of the thin d isc is a func- 
tion o f both age and the parameters zo, to, and a (iGirardi et al.l 
l2005l) . The parameters are determined by fitting the numerical 
Just-JahreiB model relationship of the thin disc. We fit the thin 
disc parameters in exponential and sech^ form to the Just-JahreiB 
model, respectively. The latter fit results in a smaller (namely 

4.7 X 10"^) than the former one (8.3 x 10"^), hence we choose 
the sech^ form as the profile of the thin disc. The derived local 
surface density of the thin disc is Sd(0) = 32.26 Mqpc"^, which 
is 10% higher than the surface density of the thin disc of the 
Just-JahreiB model. 

The halo of the Just-JahreiB model is fitted by the 
TRILEGAL form in the oblate r^^^ spheroid. The best-fit param- 
eters are local density ph(0) = 0.001576 Mqpc"^, oblateness 

= 0.65, and eff'ective radius r^^ = 1 1994 pc. 

To take into account the estimated total local density of all 
stars ever formed, the present-day densities of the Just-JahreiB 
model were corrected accordingly (see Sect. l3.1b . 



4.2.2. Star counts 

The second row in Figs. |4] and \5\ shows the star count results 
of TRILEGAL based on the Just-JahreiB parameters. The lumi- 
nosity functions are of poorer quality than that for the default 
parameter set 1, because the halo is now over-represented by 
about 0.2 dex. In the CMD, the thick disc fills the minimum in 
the centre of the Hess diagram seen in the default model. The 
F-turnoff' of the thick disc is still slightly too red (metal-rich). In 
the faint (distant) regime, there are significantly too many stars 
in all components. 

The overall fit of the Hess diagram with (x^) = 214.65 is 
worse than for the default model owing to the poor fit in the 
faint magnitude range. The median absolute relative diff'erence 
in Hess diagram Amed is 0.2993. 



Table 4. The fitting results of the thin disc sub-populations in the 
Besan^on model. 



Age range 


t 


hi 


P/(o) 


(Gyr) 


(Gyr) 


(pc) 


Mopc-3 


0-0.15 


0.075 


32.5 


0.0019 


0.15-1 


0.575 


97.8 


0.0066 


1-2 


1.5 


120.0 


0.0073 


2-3 


2.5 


195.9 


0.0043 


3-5 


4.0 


273.5 


0.0049 


5-7 


6.0 


328.3 


0.0030 


7-10 


8.5 


370.9 


0.0020 



Notes. Col. 1 lists the age range of each sub-population in the Besangon 
model. Col. 2 is the simple arithmetic mean age in each bin. We fit the 
vertical profile of each sub-population via two parameters: scale height 
hi (Col. 3) and local density p/(0) (Col. 4). The relationship of t to hi is 
determined by fitting a power-law function. 



4.3. Besangon input (set 3) 
4.3.1. Model parameters 

The thick disc of iRobin et al.l (l2003l) is described by a parabolic 
shape up to \z\ = 400 pc and an exponential profile beyond this 
height. 

We fit this piecewise-defined function using a sech^ and 
an exponential function, respectively. The exponential function 
with a smaller fitting error is adopted in Col. 5 of Table O The 
scale height of 910 pc and a local density of 0.0015 MQpc"^yield 
a local surface density X(0)=2.82 Mopc"^. 

There are seven sub-populations in the thin disc of the 
Besangon model. Each sub-population has its own age range, 
metallicity, and spatial vertical distribution. We fit scale heights 
and local densities of each sub-population (see Table |4l). 
Exponential profiles yield smaller values. 

In the next step, the scale heights as a function of age are 
fitted by the TRILEGAL law. The best-fit parameters zo, to, and 
a are 13.61 pc, 0.02016 Gyr, and 0.5524, respectively. The to- 
tal local density Pd(0) is the sum of p/(0) of each population, 
namely 0.030 Mqpc"^. The derived local surface density of the 
thin disc is 25.20 Mqpc"^. The density was corrected by the fac- 
tor 1 /geff = 1.53 to obtain the mass of ever formed stars. 

The local density of the halo in the Besangon model is de- 
termined by the SFH and the IMF. The adopted value is much 
lower than our local calibration established using SDSS data. 
We therefore fit the halo of the Besangon model with the effec- 
tive radius, oblateness, and local density of the TRILEGAL de 
Vaucouleur profile in two steps: 

1 . We fit the vertical profile to find the eff'ective radius and the 
oblateness; 

2. We use the local density as a new fitting parameter to repro- 
duce the halo star counts of the SDSS data and then use this 
value to fit a model of the parameter set 3 (Table O. 

In the first step, we find the best-fit parameters = 0.7 land 
Teff = 39.436 kpc. In the second step, we determine the local 
density of the halo by fixing and r^^ and fitting the SDSS data 
in a restricted area of the Hess diagram where the halo dom- 
inates (i.e., in the dashed triangle of Fig. O see Section 14.61) . 
The best-fit local density of the stellar halo is ph(0) = 0.00009 
Mqpc"^, which is about ten times higher than the value used in 
the Besangon model. 
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4.3.2. Star counts 

The third row in Fig. [4] and [5] shows the star count predictions of 
TRILEGAL based on the Besangon parameters. The luminos- 
ity functions appear to disagree significantly with these predic- 
tions at bright magnitudes primarily owing to the presences of 
the thin and thick discs, but be a reasonably fit at the faint end. 
In the CMD, the thin and thick disc regime are strongly under- 
represented in the model. The halo with the original local density 
also appears to be under-represented. 

The overall fit to the Hess diagram with (^) - 207.60 is 
tighter than that of parameter set 2 but worse than that of the 
default set (set 1). The Amed is 0.3445, which is worse than for 
sets 1 and 2. 

This result is not a judgement on the Besangon model com- 
pared to the Just-JahreiB or TRILEGAL model. It is a measure 
of the flexibility of the TRILEGAL model. A comparison of the 
Besan9on model with the NGP data is later discussed in more 
detail. 

4.4. Jt/r/c /npuf (se\ 4) 
AAA. Model parameters 

We do not consider the scale height as a fun ction of the age of 
the subpopulation in the Juric case {a = 0). lJuric et al~ l (l2Q08h 
fit the basic parameters of their Galactic model via star count 
observations from the S PSS. Only the fina l parameters of their 
discs are considered here. I Juric et al.l (l2008l) describe the vertical 
profiles of the thin and thick disc by simple exponentials with 
/Zd = 300 pc and pd(0) = 0.029 Mopc"^for the thin disc and 
hx = 900 pc, pt(0) = 0.0034 Mopc'^for the thick disc. This 
corresponds to 12 % of the local density of the thin disc given by 
lJuric et all (l2008l) . The derived local surface density of the thick 
disc is Et(0) = 6.12 Mqpc"^. The density has been corrected to 
the mass of ever formed stars with the factor geff for both the thin 
and thick discs. 

The halo density profile is fitted by TRILEGAL with three 
parameters. The best-fit parameters are an oblateness of = 
0.59, a local density of ph(0) = 0.000145, and an eff'ective radius 
of reff = 18304 pc. 

4.4.2. Star counts 

The bottom rows in Fig. Hand [5] show the star count results of 
TRILEGAL based on the Juric parameters. The luminosity func- 
tions are slightly more clisely reproduced than for set 3. They 
show a significant deficiency at bright magnitudes dominated by 
the thin and thick disc and too many stars at the faint end. In the 
CMD, the fit of the thick disc regime and the halo is reasonable 
in the model. At the red end, the bright part is under-represented 
and the faint end is over-represented. 

The overall fit of the Hess diagram with (x^) = 137.88 is the 
best of the four parameter sets but has the worst Amed (0.4639). 

4.5. Discussion 

An overall comparison of sets 1 to 4 shows that the default pa- 
rameter set fits best the luminosity functions in the g and r band. 
On the one hand, this confirms the fitting procedure used in 
TRILEGAL to determine these parameters mainly by compar- 
ing luminosity functions in diff'erent fields. On the other hand, 
it is remarkable that the thick disc component, which becomes 
obvious in the Hess diagram analysis, is extremely poorly fit. 



All other model adaptations fail to closely reproduce the lu- 
minosity functions. No model with a thin disc, thick disc, and 
halo can reproduce the shallow regime of the luminosity func- 
tion in the range 17mag < g < 19mag. The analysis of star 
counts in the CMD shows that disagreements between the data 
and models exceed ±50% in some regions of the Hess diagram. 
In this sense, all models fail to provide a good fit of the CMD. 

The default model shows that luminosity functions alone 
are inappropriate for determining the structure of the Milky 
Way owing to the large distance range in each component. 
Nevertheless, the total value for the Hess diagram fit is rel- 
atively small compared to the other sets. 

In the next section, we try to optimise the TRILEGAL pa- 
rameters to find the best model for the CMD at the NGR 

4.6. Optimisation (set 5) 

To find the best-fit TRILEGAL model, it is necessary to identify 
the best-fitting parameters. Because there are so many parame- 
ters that can be adjusted within TRILEGAL, we cannot deter- 
mine them simultaneously. However since the thin disc, thick 
disc, and the halo do not completely overlap in the (g - r,g) 
CMD, it is possible to optimise the parameters for the halo, thick 
disc, and thin disc sequentially. 

Fig- gives the schematic diagram for distinguishing areas 
where only the halo or both the halo and thick disc contribute 
significantly to the star counts, and where only the thick and thin 
discs contribute in this manner. We define the separating lines in 
terms of the relative densities of these three components based 
on the Just-JahreiB model. In the resulting triangle (dashed line), 
the halo dominates over the other two components because the 
relative density of the discs to the halo is less than 1%. In the 
area below the solid line, the thin disc can be neglected, because 
the relative density of thin disc to thick disc is less than 1%. We 
fit first the halo parameters in the triangle. We then fix the halo 
parameters and fit the thick disc parameters below the full line by 
including the star counts of the halo. Finally, we investigate the 
thin disc parameters by taking into account all three components 
in the full CMD. 

4.6.1. Halo fitting 

Locating an area with halo stars but very few disc stars in the 
CMD is the first step. Our choice is marked by the dashed tri- 
angle in Fig. [6l The discs stars in this selected region can be 
neglected, i.e., (pthin + Pthick)/Phaio < 0-01 as estimated by the 
Just-JahreiB model. The boundary conditions to determine the 
halo parameters are g - r > 0.2 and g > 15.5 -h I2.5(g - r). 
We choose an oblate r^^^ spheroid and search for a best fit in the 
parameter ranges 

0.00005 Mo/pc^ < Qh(0) < 0.00035Mo/pc^ 
2000 pc < reff < 30000 pc, 
0.57 <qh< 0.70. 

The fitting result is shown in Fig. [71 The best-fit halo pa- 
rameters are listed in the last column (set 5) of Table [3] with 
(X^) = 4.86. 

4.6.2. Thick disc fitting 

We fix the parameters of the halo determined in Sect. 14.6.11 to 
fit the parameters of the thick disc in the area dominated by 
the thick disc and halo (below the solid line determined by 



11 



S. Gao et al: Comparison of the Milky Way model and SDSS data at the NGP 



Thin + thick discs dominate 





Fig. 6. A schematic diagram of the CMD illustrating the distinc- 
tion between thin and thick disc and halo stars. The dashed trian- 
gle represents the area where the halo dominates. The solid line 
is the upper boundary below which the thin disc contribution is 
negligible. 



g > 16.5 -h 3.64(g - r) in Fig. O. In this area, only relatively 
few thin-disc stars exist, i.e., Pthin/(Pthick + Phaio) < 0.01. 

The sech^O and exp() prescriptions are tried to fit the vertical 
profile of thick disc with the parameters scale height ht and local 
calibration pt(0) within the following ranges: 

0.0010 Mo/pc^ < Pt(0) < 0.00035Mo/pc^; 
400 pc < ht < 900 pc. 

The minimum (x^) = 16.22 was found for the sech^O profile 
with the best-fitting parameters listed in Col. 7 of Table [S] The 
minimum (;^2)ofthe fit using an exponential function is 0.312. 

The fit result, shown in Fig. [8l shows on the one hand the 
homogeneous quality of the fit, but on the other hand a system- 
atic trend along the MS of the thick disc with too few stars in the 
G-dwarf regime and too many K and M dwarfs. 

4.6.3. Thin disc choice 

For the thin disc fitting, we now fix the halo and thick disc as 
determined above. For the thin disc, we cannot get a best fit for 
all the free parameters. On the basis of the analysis of the advan- 
tages and disadvantages of the thin disc models of for the param- 
eter sets 1, 2, 3, and 4, we adjust the parameters iteratively to get 
a better result. All iterations to improve the parameters in order 
to find an optimal Milky Way model must be done by hand. The 
power-law index a in the scale-height function changes the frac- 
tion of diff'erent sub-populations as a function of height above 
the mid-plane. In the Hess diagram, this aff'ects the balance of 
densities between g-r > I (faint, nearby) and 0.2 < g - r < 0.4 
(bright, more distant). The values for zo and pd,o aff'ect the to- 
tal density of the thin disc. These are the general directions for 
adjusting the parameters. The value of to was not changed. 

We performed these investigations for both vertical profile 
shapes. The sech^ profile results in a better fit with a minimum 
(X^) = 1.33. (The exponential profile leads to a best fit with 
<X^} = 2.95). 
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Fig. 7. Hess diagram, relative diff'erences, and x'^ distribution 
(top to bottom) of the best-fit halo in TRILEGAL for parame- 
ter set 5. The upper panel is the Hess diagram of the halo of set 
5. The white thick lines enclose the triangle mentioned in Fig. (6] 
Two lower panels are the relative diff'erence and x^ distribution 
restricted in the halo area. The coding is the same as in the rows 
of Fig. [5] (A colour version is available on-line.) 



4.6.4. Model properties 

The Hess diagram analysis for input parameter set 5 is plotted 
in Fig.|9l The;^^ distribution is smoother than for the other four 
sets, but the systematic features in the relative deviations of data 
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Fig. 8. Same as in Fig. [7] but for thick disc plus halo stars and 
the parameter set 5. The upper panel is the Hess diagram of the 
halo plus the thick disc of Set 5. The white thick line in Fig. [6l 
shows the area in which the thin disc is negligible. The two lower 
panels are the relative difference and;^^ distribution restricted in 
the thick disc area. (A colour version is available on-line.) 



from model do not vanish. These are as follows(see middle panel 
of Fig. O: 

1. The dark vertical stripe at g - r = 0.2 is a sign of too 
many blue thin-disc stars (A stars), which do not show up 
in Figures [8] and [TJ The constant SFR overestimates the frac- 
tion of young stars. 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

g-r (mag) 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 
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Fig. 9. Hess diagrams (top), relative differences (middle) be- 
tween the Hess diagrams of the parameter set 5 and the SDSS 
observations, and the;^^ diagrams (bottom). The notation is the 
same as in Fig.[5l (A colour version is available on-line.) 



2. The white island at g - r = 0.35 and 15 < g < 17 mag points 
to a mismatch of the thick disc turn-off stars. There are too 
few nearby stars owing to the flat sech^ profile and the turn- 
off is too red. 

3. The faint red end with a colour-magnitude range of g - r > 
1.0 and g > 18 mag shows too many M dwarfs. The reason 
may be that the g-r colour of the thin disc M dwarfs is 
slightly too blue or that the IMF allows too many stars at the 
low mass end. 



13 



S. Gao et al: Comparison of the Milky Way model and SDSS data at the NGP 





3 


1 




















CD 


1.5 


■o 












1 


-z. 


o 




oT 


0.5 













16 18 

g (mag) 




14 



16 18 

r (mag) 



20 



Fig. 10. The magnitude distributions of the SDSS observations 
and the TRILEGAL simulations generated with parameter set 
5 in the g (top) and r (bottom) band. The x- and y-axes are 
similar to the bottom panel of Fig. [T] The error bars (blue) are 
from the SDSS DR7 data with a colour limit of 0.2 < (g - r) < 
1.2. The black solid line represents the model star counts from 
TRILEGAL set 5. The dotted (red), dashed (cyan), and dot- 
dashed (green) line stand for the thin disc, thick disc, and halo 
of TRILEGAL set 5, respectively. (A colour version is available 
on-line.) 



4. The white islands at the bright red end are another hint that 
the red giants of the thick disc and halo are too blue (see also 
Figs.SandE]). 

The magnitude distribution in the g and r bands (see Fig. 
[Tot shows larger discrepancies than for the default Set 1 (top 
row of Fig. 131). The shoulder at 17mag < g < 19mag is still 
missing. This is again a hint that the Hess diagram analysis is 
a much more powerful constraint of the stellar population than 
luminosity function fitting. 

To quantify the impact of giants on the Hess diagram, we 
select the giants using the surface gravity logg < 3.8. The abso- 
lute and relative contributions of giants simulated by parameter 
set 5 are shown in Figure [TTJ The giants of ~ 30% stars exist 
in the Hess diagram. The rigid line of giant counts is close to 
g-r = 0.5~0.7. The giants in the thin disc have a higher rela- 
tive contribution, though the absolute counts are too low for the 
thin disc. 

The parameters of the optimized model are listed in Table 
[3] (set 5). Compared to the default model (Set 1), the final Set 
5 has a lower local surface density, a smaller maximum scale 
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Fig. 11. The fraction of the contribution by modelling giants 
of parameter set 5 in the TRILEGAL model. The upper panel 
shows the giant counts simulated by set 5 of the TRILEGAL 
model in the NGP field. The colour code covers the star number 
density from 1 to 10 per square degree, 1 mag in g, and 0.1 mag 
mg- r. Lower panel: The relative fraction of giants to all stars 
in the model. The colour code is from 0% to 40%. (A colour 
version is available on-line.) 



height /Zz,max of the thin disc, and a shallower nearby halo profile 
of lower local density, because the old, extended thin disc and 
the halo do not need to account also for the thick disc stars. 

Compared to the Just-JahreiB model with i^x^) - 4.31 
(I Just et al.l I2OIII) . the fit of Set 5 is much worse. This is also 
obvious from the median relative deviations Amed — 

0.2576. The 

mismatch of the isochrones in the ugriz filters alone cannot ac- 
count for these deviations. The simple structural forms of the 
thin and thick discs are the main reasons for the large discrepan- 
cies. 



5. Simulations and results of the Besangon model 

The Besangon model allows no variations in the choice of free 
parameters. We performed a multi-colour analysis of the predic- 
tions in diff'erent colours, as given in Table [5] We selected all 
stars in the NGP field from DR7 within the colour and apparent 
magnitude ranges listed in Table [5|3 and with magnitude errors 



^ The mock catalogues of the thin and thick discs and the halo and 
of their sum are generated by Monte Carlo procedures in different pro- 
cedures. The number of stars in the sum of the three components are 



14 



S. Gao et al: Comparison of the Milky Way model and SDSS data at the NGP 
Table 5. Properties of the observed and simulated CMDs. 



dia^ 


^ram 




colour 


magnitude 


Nobs 


N mod 


^ {Nthin 


+N thick 


+Nhalo) 


N modi N obs 


<x' > 


^med 


(u - 


8^ g) 


[- 


-0.2,2.0) 


[14.0, 20.5) 


238 724 


261 205 


43 285 


133 603 


84550 


1.09 


555.43 


0.5256 


(g- 


r,g) 


[- 


-0.2,1.2) 


[14.0, 20.5) 


279 174 


312 093 


94 590 


134 209 


83 587 


1.12 


408.15 


0.2724 


(r- 


Ur) 


[- 


-0.2, 1.2) 


[14.0, 20.5) 


479 913 


433 533 


145 784 


183 777 


103 955 


0.90 


222.62 


0.2595 


(i- 


zj) 


[- 


-0.2,0.5) 


[14.0,18.5) 


220 144 


237 249 


97 560 


114468 


25 236 


1.08 


160.60 


0.2029 



Notes. Col. 1 gives the colour and magnitude of the CMD. Cols. 2 and 3 list the ranges in colour and apparent magnitude, and Cols. 4 and 5 show 
the total number of stars in the data {Nobs) and model {N^od) for each CMD. The stellar numbers of the three components are listed as Cols. 6 to 8 
{Nthin, Nthick and Nhaio)- Col. 9 shows the ratio of star numbers of the model to the observations. Col. 10 lists the < > of the Hess diagrams. The 
last column lists the median absolute relative difference t^rmmed of each Hess diagram for the region redder than the turn-off. 



smaller than 0.2 mag in each filt er. The selection of the field and 
objects is similar to that used in I Just et aP (1201 ll) and Section 4. 

5.t. Isochrones 

Since the thick disc and the halo are modelled in the Besan^on 
model by simple stellar populations with a single age and a rel- 
atively narrow metallicity range (11 Gyr and [Fe/H] = -0.78 ± 
0.30 for the thick disc, 14 Gyr and [Fe/H] = -1.78 ± 0.50 for 
the halo), we can compare the isochrones with the fiducial se- 
quen ces of star cluster s with corresponding properties (see Table 
[T|. In I An et all (l2009b . a detailed analysis of isochrone fitting to 
fiducial s equences i n the u ^riz filters of globular and open clus- 
ters from lAn et al.l (l2008l) is presented. The authors included a 
comparison with c l usters o bserved with the CFHT-MegaCam by 
IClem etal.1 (l2QQ8h . lAn eTa l. (2009) found that the clusters show 
a systematic offset relative to the Yale Rotating Evolutionary 
Code (YREC) and MARCS models owing to a zero-point prob- 
lem. We use five old globular clusters to cover the metallicity 
range of the thick disc and halo. 

We plot CMDs using the absolute magnitudes Mg, Mr, and 
Mi from -3 mag to 15 mag and compare the CMDs to the fidu- 
cial isochrones of the globular clusters presented in Table [T] 
There is a considerable offset of 0.3 mag m(u- g) between the 
fiducial isochrones and the simulated CMDs (top row of Fig.[T2|, 
which is a hint of a zero-point problem in the u band. This large 
offset explains why the Hess diagram of the u-g completely fail 
to the actual data (see below). 

In Fig. [121 the left panels show the comparison for the thick 
disc, whereas the right-hand panels are for the halo. In the dif- 
ferent CMDs, the fiducial isochrones of the globular clusters are 
overplotted: M5 and M71 are used to represent the thick disc; 
and M5, NGC 4147 and M92 are used for the halo (plotted in 
sequence from the turn-off point from the left to right in the right 
panels). 

In the other CMDs of the thick disc, there is a smaller, 
but still significant, offset in the blue colour for the MS in the 
Besangon model compared to the fiducial isochrones. The metal- 
licity of the clus ter M7 1 is closest to the metallicity of the thick 
disc adopted in iRobin et al.] (l2003h . The matches of M5 with 
the lower metallicity are obviously better than for M7 1 , though 
there is still an offset(~ 0.2 mag in the u - g) between the fidu- 
cial isochrones of M5 and the sequence of the thick disc of the 
Besangon model. The isochrones of the thick disc correspond to 
fiducial isochrones with a metallicity below -1.0. 

For the halo, the offsets are similar but slightly larger. The 
isochrone with [Fe/H] = -2.28 is closest to the model but cor- 



therefore not exactly the same as the number of stars in the general 
simulation. 



responds to the lower boundary of the metallicity range in the 
Besangon model. 

5.2. Luminosity functions 

In Fig. [131 the luminosity functions of the SDSS data (with er- 
ror bars) and the simulated data are compared in all five filters. 
The full (black) step-function shows all stars of the Besan^on 
model with contributions from the thin disc (dotted red lines), 
thick disc (dashed cyan lines), and halo (dot-dashed green lines). 
The model fits the data well over a wide magnitude range in all 
filters. A few systematic features can be seen: 

- In the u, g, /, and z bands, the predicted number densities are 

slightly higher than in the actual data. At the brighter mag- 
nitudes (g < 15 mag), these differences are significant in all 
five filters. The thin disc dominates at the bright end where 
the differences between the predictions and the data show 
the largest deviations (0.2 dex to 0.5 dex). These stars corre- 
spond to F and early G stars (Mg = 4 to 5 mag) of the MS o f 
the thin disc with distances above ~ 500 pc (see Sect. 15.3b . 
The reason for the discrepancy may be a higher mass density 
of the thin disc at large z or a higher fraction of early-type 
stars caused by the IMF or the SFH. 

- Only in the r band, the star counts are underestimated at r > 

- The thin disc dominates only at the bright end, but contributes 

significantly up to the faint magnitude limit in the r, /, and z 
bands. 

- The contribution of the thick disc dominates over wide mag- 

nitude ranges showing that the thick disc model parameters 
are very important in attempting to reproduce the luminosity 
functions. 

- The halo contributes significantly only at the faint end of the 

u, g, and r band diagrams. 

5.3. Hess diagrams 

We calculate the Hess diagrams of the SDSS data and the 
Besangon model in the same way. The star counts are binned 
in each CMD in intervals of 0.05 mag in colour and 0.5 mag 
in magnitude and smoothed in steps of 0.01 mag. Numerical 
counts are given on a logarithmic scale covering a range of 1 
to 100 per deg^, 0.1 mag in colour, and 1 mag in luminosity. 
The relative differences between the model and the data (i.e., 
(data - model)/model) are smoothed only in luminosity and 
cover a linear range of -1.0 to +1.0. We investigate star counts 
in the four CMDs (u - g, g), (g - r, g), (r - /, r), and (/ - z, 0- To 
show and compare each stellar population in all the CMDs, the 
colour and magnitude ranges are different in the different CMDs 
(see columns 2 and 3 of Table [5l). 
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Fig. 12. The colour-absolute magnitude diagrams of the thick disc and halo of the Besangon model and the fiducial isochrones of 
globular clusters in {g - r, Mg). The Hess diagrams show the stellar distributions in {u - g, Mg), (g - r, Mg), (r - /, Mr), and (/ - z, Mi) 
from top to bottom for t he thick d i sc (le ft column) and halo (right column) simulated by the Besangon model. The overplotted 
thick lines (white) are the I An et aL I (l2008l) 's fiducial isochrones based on SDSS observations of globular clusters. M5 and M71 are 
overplotted on the thick disc panels. M5, NGC 4147, and M92 are overplotted on the halo panels. (A colour version is available 
on-line.) 
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Fig. 13. The distributions of apparent magnitude of the Besan9on model in the r, /, and z bands around the NGP. The thin disc, 
thick disc, and halo are plotted as dotted lines (red), dashed lines (cyan), and dot-dashed lines (green), respectively. The solid lines 
(black) represent their sums. The error bars (blue) show data from SDSS DR7 for comparison. The samples of the u and z band data 
are selected in colour intervals of -0.2 <(u-g)< 2.0 and -0.2 < (/ - z) < 0.5, respectively. The other colour ranges extend from 
-0.2 to 1.2 (see second column in TableQ. (A colour version is available on-line.) 



In Figures [141 -[T71 the results are shown for all CMDs. The 
upper panels show the Hess diagrams of the data, the model, and 
their relative differences (from left to right). In the lower panels, 
the contributions of the thin disc, the thick disc, and the halo 
(from left to right) as given by the Besangon model are shown 
with the same colour coding. 



We first discuss t he results in (;?-r, ^) f or a direct comparison 
to the Just-JahreiB (|JusL&MirdB|[2^ and the TRILEGAL 
(iGirardi et al .120051) models. The top right panel in Fig. [T5l shows 
that the relative diff'erences for most regions are smaller than 
50%. There are four regions in the CMD exceeding that level: 



- The deep minimum at the lower right corner, where the star 

count predictions exceed the observations by more than a 
factor of three (compare the middle and left panels in the top 
row) is due to M dwarfs in the thin disc model (lower left 
panel). It seems that the (g - r) colour of the M dwarfs is too 
blue. 

- The F turn-off' regions of the thick disc (lower middle panel) 

and the halo (lower right panel) are also blue- shifted by ~ 
0. 1 mag in {g - r), leading to the negative-positive features at 
g = 16 and g = 20 mag, respectively, in the diff'erence plot. 

- At the bright end (g = 14 to 15 mag), the star counts are deter- 

mined by MS stars of the nearby thin disc, as well as bright F 
turn-oflf stars and giants in the thick disc. The thick-disc gi- 
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Fig. 14. Hess diagrams of the SDSS data and the Besan9on model in (u - g, g). The upper three plots show the Hess diagram of the 
SDSS data, the model, and the relative difference between the data and model (from left to right). The lower panels show the Hess 
diagrams of the three simulated components of the Besan^on model: thin disc, thick disc, and stellar halo (from left to right). All 
Hess diagrams cover a range of stellar densities from 1 to 100 per square degree per magnitude in luminosity and 0.1 magnitude 
in colour (purple to red in the electronic version; black means yet lower number densities). The difference plot (top right) covers a 
range of relative differences from -100% to +100% (purple to red in the electronic version; black and white mean larger negative 
or positive differences, respectively). The ranges in apparent magnitude and colour of these six plots are listed in columns 2 and 3 
of Table [3) . (A colour version is available on-line.) 
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Fig. 15. Same as Fig.[T4]but in (g - r, g). 



ants are responsible for the excess of stars on the blue side of 
the diagram. The maximum at the top right corner is a sign 
of missing thin-disc K dwarfs or red giants in the thick disc. 
The faint blue part at g > 19 mag and ^ - r = 0.2 to 0.4 mag 
is dominated by the halo. The density of the halo is strongly 
under-represented by the model. This cannot be seen in the 
luminosity function (see Fig. [13]), since it is compensated 



for by the over-estimated number of thin disc M dwarfs at 
g-r < 1.2. 

The Besangon model confirms that the faint, very blue plume 
at g - r < 0.2 cannot be reproduced by stellar populations 
of the Milky Way and is most probably dominate d by mis- 
identified extragalactic sources (ILupton et al .120011) (see also 
iSchultheis et al...2006i) . 
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Fig. 16. Same as Fig.[T4]but in (r - /, r). 
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Fig. 17. Same as Fig.[T4]but in (/ - z, /)• 



In the Just-JahreiB model, the median relative deviation in 
the Hess diagram is Amed = 5.63%. This corresponds to a re- 
duced value, measured by the star counts in logarithmic scale, 
of 4.31. In the Besan^on model, typical deviations are five times 
larger with a corresponding;^^ = 408.15. The median relative 
deviations of the diff'erent colour diagrams are between 20% 
and 53%. In the CMD (Hess diagram), the star count distribu- 
tion at a given colour is mainly determined by the geometrical 
structure of the components, whereas the distribution at a given 
magnitude is strongly influenced by the colours of the underly- 
ing stellar populations. Since all strong features in the relative 
difl'erences between the data and model in Fig. [15] show a hor- 
izontal structure, the main reason for this disagreement with in 
the Besangon model are the mis-matching properties of the syn- 



thetic stellar populations (see Fig. [12]). In the Just-JahreiB model, 
the MS properties of the stellar components fit by construction 
in the solar neighbourhood. 

In the TRILEGAL model (Section 4), the colour off'sets in 
the synthetic stellar populations are much smaller and do not 
completely dominate the deviations of the data from the model. 
For features in the Hess diagram other than the ones just dis- 
cussed, the absolute relative diff'erence has a median value of 
0.2576 corresponding to = 104.19. 

In case of the Besangon model, the overall diff'erences be- 
tween the data and model arising from deviations of the density 
profiles are large and cannot be quantified here. 

In (u-g, g) (see Fig. [14]) diagram the diff'erences between the 
data and the Besan^on model are much larger than in the (g-r, g) 
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plot. The main reason is that the isochrones are even more blue- 
shifted by ~ 0.3 mag in (u - g) (see Fig. [121). In addition, the 
features in the number density distributions of the Hess diagrams 
are much wider in colour in the observations than in the model. 

In the (r - /, r) (see Fig. [161) diagram, the thin disc dominates 
for r - / > 0.8, while the thick disc dominates for 0.3 < r - i < 
0.7. That the largest significant diff'erence between the model 
and the observations is interestingly located at the interface of 
the populations of the thin and thick disc is caused by a gap at 
0.75 < r - i < 0.85 and 20 < r < 20.5. Another feature is 
a pronounced gap at r - / ^ 0.45 in the MS of the thin disc 
(lower left panel), which is much weaker in the other colours. 
These minima correspond to the relative number counts along 
the MS of the thin disc. At red colours (low-mass end), the star 
numbers are determined by the IMF, whereas in the blue (high- 
mass end) there is a strong impact of the SFH owing to the finite 
MS lifetime of the stars. 

The Hess diagram of (/ - z, /) (Fig.[T7l) represents the best fit 
of the Besangon model and SDSS data. In most regions of the 
CMD the median absolute diff'erence is about 0.20. 

Altogetherour multi-colour analysis confirms that the syn- 
thetic population properties first of all need to be improved in 
the Besangon model. In a second step, the IMF, SFH, and the 
spatial structure should be investigated with higher precision. 

6. Summary 

We have compared star-count predictions generated by the 
TRILEGAL and Besangon models with the observed sample in 
the NGP field (b > 80°) from the SDSS DR7 photometric stellar 
catalogu e using de-reddened rn agnitudes adopting the extinction 
values of ISchlegel et al.l (Il998l) . We measured the quality of the 
agreement between diff'erent models and the data by perform- 
ing 8iX^ determination of the Hess diagrams, i.e., the star count 
distribution in the CMD. As a comparison, we referred to the 
Just-J ahreiB model (model A in I Just & Jahreii3l[2010l: I Just et al.l 
I2OIII) with (x^) = 431 corresponding to a median relative devi- 
ation between the data and model of Amed = 5.63% in (g - r, g). 
Additionally, fiducial isochr ones of five glob ular clusters, based 
on the SDSS observations o f lAnetalJ (120081) combined with the 
cluster parameters of iHarrisI ( 1996h . were used to test the syn- 
thetic stellar populations of the thick disc and halo modelled by 
the TRILEGAL and Besan^on models. 



6.1. The TRILEGAL model 

We studied with TRILEGAL five diff'erent models simulating 
the NGP region and compared their output to observational pho- 
tometric data of this region from the SDSS. The web inter- 
face of TRILEGAL allowed us a restricted choice of input pa- 
rameters and functions for the Galactic components. We com- 
pared the default input set (lacking a thick disc component) 
with best-fit adjustments for three other Galax y models from 
the li terature, namely the Just-JahreiB mod el (jJust & JahreiBi 
201 oh. the Besang on model (iRobin et al.l2003h . and the model of 
Juric et all (l2008l) . We found that the default set reproduces best 
the luminosity functions in the g and r bands. The (y^) = 151.98 
from the Hess diagram analysis is slightly worse than that based 
on the Juric model with (x^) = 137.88. All these models show a 
typical discrepancy Amed of 26% to 46% in the Hess diagram, 
which is significantly larger than the disagreements from the 
Just-JahreiB model. 

A comparison of the colour- absolute magnitude diagrams of 
the thick disc and halo with the observed fiducial isochrones 



of globular clusters (lAn et al ] l2008h shows that the MSs and 
the F-turnoff' regions of the stellar evolution library are roughly 
consistent in the g and r filters. The small colour off'set of 
g - r ~ 0.05 mag may be due to age diff'erences among the 
components. However, the giant branches show large system- 
atic disagreements. The red colour of the K giants observed in 
the globular clusters are not reproduced. 

Since the parameter space was too large to find the absolute 
minimum (x^) with TRILEGAL, we determined the parameters 
for halo, thick disc, and thin disc sequentially by fitting smaller 
areas of the CMD. Finally we found an optimised set 5 with 
(X^) = 1.33 and a median deviation Amed of 26%, which is five 
times larger than in the Just-JahreiB model. The parameters of 
this model are listed in the third-last column of Table [3l The 
reason for the larger deviations is mainly the restricted choice 
of thin and thick disc parameters (for a constant SFR and the 
adopted vertical density-profile shapes). 

The TRILEGAL code is a powerful tool for testing Milky 
Way models against star count data. The investigation of star 
number densities in the CMD (Hess diagrams) are particularly 
helpful in distinguishing the impacts of the diff'erent stellar pop- 
ulations. 

Some more progress may be possible by varying the thin disc 
parameters further, but the main disagreements are due to funda- 
mental limitations of the models. In the TRILEGAL model, the 
range of input parameters and functions is quite restricted. More 
freedom in choosing the SFH, the vertical profiles of the thin and 
thick discs and for halo profile would be desirable when simu- 
lating high Galactic-latitude data. The comparison with fiducial 
star cluster isochrones has additionally shown that an improve- 
ment in the predicted colours of stars for late stages of stellar 
evolution in the ugriz filter system is desirable. 

6.2. The Besangon model 

We investigated in a multi-colour analysis the properties of the 
Besangon model by comparing with SDSS data of the NGP field. 
We generated mock stellar catalogues with the Besangon model. 
The output in u*, g\ r\ i\ and photometries was converted 
into the SD SS u, ^, r, i, and z s yste m using the linear rela tions 
provided bv lTucker et"aD (l2006h and lJust & JahreiBi (l2008h . 

The comparison of the colour- absolute magnitude CMDs of 
the model output of thick disc and halo with the fiducial globu- 
lar clusters covering the respective metallicity ranges, uncovered 
significant diff'erences. The MS and turnoff' colours are blue- 
shifted by from 0.05 mag to 0.3 mag, increasing from the red to 
the blue colour bands. The large off'set inu- gis mainly caused 
by a zero-p oint problem in the u' filter of the CFHT-MegaCam 
reported by I An et aP (l2008l) . The giant branches are fitted well 
but the models do not extend far enough to the red end of the 
CMDs. 

A comparison of the luminosity functions in the diff'erent fil- 
ters shows that there is a good general agreement between the 
data and the model. Only for magnitudes < 15 mag, the star 
counts are overestimated. In the r band, there is a deficit of stars 
at the faint end at r > 19 mag. The contributions of the diff'er- 
ent components (thin disc, thick disc, and halo) as a function 
of apparent magnitude were plotted to identify the dominating 
populations. The fraction of giants are generally of the order of 
a few percent, which may be neglected for the sake of our tests. 

The Besangon model output and the diff'erences from the 
data were plotted in Hess diagrams of four colours [(u - g, g), 
(g -r,g),(r- /, r), and (/ - z, i)] . To analyse the impact of the dif- 
ferent Galactic components, the individual Hess diagrams of the 
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three components were also shown. The relative difference be- 
tween the data and the model cover a range larger than - 100% to 
100% in each Hess diagram. Except for some well-defined fea- 
tures caused by the colour off'set of the isochrones, the median 
relative diff'erences in the Hess diagrams are between 0.2 and 0.5 
(ignoring the strongly perturbed (u - g, g) diagram). The largest 
discrepancies in the Hess diagrams are the follows: 

- The MS turn-off' regions of thick disc and halo are underesti- 

mated in terms of the observed colours, meanwhile they are 
overestimated in terms of the modelled colours owing to the 
blue- shifted MS in the Besangon model. 

- The halo is strongly underestimated, which is compensated for 

in the luminosity functions by M dwarfs of the thin disc, 
which fall in the model in the investigated colour range but 
do not do so in the data. 

- The excess of observed stars in the transition between the thin 

disc and the thick disc compared to the model in terms of 
blue colours may also be due to the shifts in the synthetic 
colours produced by. 

- The prominent gap in the MS of the thin disc at r ~ 0.45, 

which is not present in the SDSS data. This feature appears 
in all of the thin disc realizations of the Hess diagrams, but 
the gap in the r - i diagram is the most significant one. 

6.3. General conclusions 

We have investigated the ability of the TRILEGAL model to re- 
produce star counts in the SDSS filters (g - r, g) and the qual- 
ity of the Besan^on model in the full-filter system ugriz for the 
NGP field. We have analysed the number distributions of stars in 
the CMDs (via Hess diagrams) instead of only us ing luminosity 
funct i ons. We have us ed the Just-JahreiB model (lJust & JahreiBi 
l2010l:lJuste'taDl201lh to quantify the quality of the fits. In com- 
parison to the Just-JahreiB model with a median deviation be- 
tween data and model of 5.6%, the TRILEGAL model with devi- 
ations of the order of 26% and the Besangon model with median 
deviations reaching from 20% to 53% show significant deficits. 
In the TRILEGAL model, the main reason for the larger devia- 
tions is the restricted choice of the density profiles for the thin 
disc, thick disc, and halo. In the Besangon model, the synthetic 
colours in the ugriz filters of the stellar populations need first of 
all to be improved, before the other properties of the Milky Way 
components can be investigated further. 

In general. Milky Way models based on synthetic stellar 
populations such as the TRILEGAL and Besangon models are 
powerful tools for producing mock catalogues of stellar popula- 
tions for comparison with data. Both models, TRILEGAL and 
Besan^on, reproduce the observed SDSS luminosity functions 
very well. However, the luminosity functions are insensitive to 
the assumed structural properties of the Milky Way components. 
To improve the modelled Hess diagrams, the synthetic colours 
of the stellar populations in ugriz filters need to be upgraded in 
both the TRILEGAL and the Besan^on models. We are looking 
forward to applying the TRI LEGAL, the Besangon, or the new 
Galaxia (ISharma et al.ll201lh tool to the full SDSS data set. 
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